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Abstract. 

Aims. A census of molecular hydrogen flows across the entire Orion A Giant Molecular Cloud is sought. With this paper we aim 
to associate each flow with its progenitor and associated molecular core, so that the characteristics of the outflows and outflow 
sources can be established. 

Methods. We present wide-field near-infrared images of Orion A, obtained with the Wide Field Camera, WFCAM, on the United 
Kingdom Infrared Telescope. Broad-band K and narrow-band H 2 1-0S(1) images of a contiguous ~8 square degree region are 
compared to mid-IR photometry from the Spitzer Space Telescope and (sub)millimetre dust-continuum maps obtained with 
the MAMBO and SCUBA bolometer arrays. Using previously-published H 2 images, we also measure proper motions for H 2 
features in 33 outflows, and use these data to help associate flows with existing sources and/or dust cores. 
Results. Together these data give a detailed picture of dynamical star formation across this extensive region. We increase the 
number of known H 2 outflows to 116. A total of 111 H 2 flows were observed with Spitzer; outflow sources are identified for 
72 of them (12 more H 2 flows have tentative progenitors). The MAMBO 1200 /jm maps cover 97 H 2 flows; 57 of them (59%) 
are associated with Spitzer sources and either dust cores or extended 1200 fim emission. The H 2 jets are widely distributed and 
randomly orientated; the jets do not appear to be orthogonal to large-scale filaments or even to the small-scale cores associated 
with the outflow sources (at least when traced with the 11" resolution of the 1200 /mi MAMBO observations). Moreover, H 2 
jet lengths (L) and opening angles (6) are not obviously correlated with indicators of outflow source age - source spectral 
index, a (measured from mid-IR photometry), or (sub)millimetre core flux. It seems clear that excitation requirements limit the 
usefulness of H 2 as a tracer of L and 6 (though jet position angles are well defined). 

Conclusions. We demonstrate that H 2 jet sources are predominantly protostellar sources with flat or positive near-to-mid-IR 
spectral indices, rather than disk-excess (or T Tauri) stars. Most protostars associated with molecular cores drive H 2 outflows. 
However, not all molecular cores are associated with protostars or H 2 jets. On statistical grounds, the H 2 jet phase may be 
marginally shorter than the protostellar phase, though must be considerably (by an order of magnitude) shorter than the prestellar 
phase. In terms of range and mean value of a, H 2 jet sources are indistinguishable from protostars. The spread in a observed 
for both protostars and H 2 outflow sources is probably a function of inclination angle as much as source age. The few true 
protostars without H 2 jets are almost certainly more evolved than their H 2 -jet-driving counterparts, although these later stages 
of protostellar evolution (as the source transitions to being a "disk excess" source) must be very brief, since a large fraction of 
protostars do drive H 2 flows. We also find that protostars that power molecular outflows are no more (nor no less) clustered 
than protostars that do not. This suggests that the H 2 emission regions in jets and outflows from young stars weaken and fade 
very quickly, before the source evolves from protostar to pre-main-sequence star, and on time-scales much shorter than those 
associated with the T Tauri phase, the Herbig-Haro jet phase, and the dispersal of young stellar objects. 

Key words. Stars: circumstellar matter - Infrared: ISM - ISM: jets and outflows - ISM: Herbig-Haro objects - Stars: mass-loss 



1. Introduction 

The southern part of the Orion constellation encompasses the 
Orion A and B Giant Molecular Clouds (GMCs), numerous 
compact and intermediate-sized molecular cores, low and high 
mass young stars (including a massive star/young OB clus- 
ter), stars at varying evolutionary stages, and dozens of op- 
tical Herbig-Haro (HH) objects and molecular outflows (see 
Peterson & Megeath| (2008) and lAllen & Davisl (2008) for re- 



views). In molecular line maps the complex extends roughly 
northwest-southeast, parallel with the Galactic plane, over 
about 13° dKutner et al.l 1977: ISakamoto et all 1994). Overall, 
the Orion A and B GMCs represent one of the richest star form- 
ing regions known. 

In extensive regions like Orion, jets and outflows can be 
used as sign-posts of on-going star formation. An abundance 
of jets points to active accretion and a young stellar popula- 
tion; a paucity of molecular outflows, in a region where near- 
and mid-IR photometry still indicate a sizable population of 
sources with reddening and excess, suggests a more evolved 
region with a fraction of pre-main-sequence stars (T Tauri and 
Herbig Ae/Be stars) observed "edge-on" through their circum- 
stellar disks. 

By using H2 2.122/im emission as a tracer of jets and out- 
flows, wide-field narrow-band images may be used to pin-point 
the locations of the youngest sources. Moreover, they allow 
one to take a statistical approach when considering questions 
about the distribution of Class 0/1 protostars and Class II/III 
Young Stellar Objects (YSOs), the interaction of these forming 
stars with their surrounding environment, the overall star for- 
mation efficiency, and the evolution of the region as a whole 
(e.g. IStankel 200 0; |Stanke, McCaughrean & Zinnecker| 2002; 
IDavis et al.l 2007: IKumar et al.l 2007 : IDavis et al.l 200"8l 

With the wide-field near-IR camera WFCAM at the United 
Kingdom Infrared Telescope (UKIRT) we have secured ho- 
mogeneously deep, sub-arcsecond-resolution H2 2.122/im and 
broad-band K images over most of the Orion A GMC. Our 
~8 square degree mosaics encompass the the molecular clouds 
known as OMC 2 and OMC 3, the Orion Nebula Cluster (ONC 
- also known as M 42), and an abundance of well-known 
low mass star forming cores, HH objects and bipolar out- 
flows spread throughout the Lynds dark cloud L 1641 (HH 1/2, 
HH 33/40, HH 34, HH 38/43, etc.). 

Our goals with these observations were to: (1) extend the 
maps of Stanke et al. (2002; hereafter Sta02) to give more com- 
plete coverage of Orion A, to search for new H2 flows, and to 
better trace the true extent of known H2 flows by searching 
for emission off the main molecular ridge; (2) provide a sec- 
ond epoch for proper motion studies, and (3) compare the ob- 
servations with extensive (sub)millimetre and mid-IR Spitzer 
observations, so that H2 jets and outflows could be associated 
with molecular cores and/or embedded protostars, the latter be- 
ing identified from mid-IR photometry. We utilise the 850 fim 



SCUBA data published by|Nutter & Ward-Thompson (2007) 



the more extensive 1200/vm MAMBO observations of Stanke 
et al. (in prep.), and the recent Spitzer observations of Megeath 
et al. (in prep.). 



This paper is structured as follows: in Sect. 3 we briefly dis- 
cuss wide-field images of regions of interest spread throughout 
Orion A and give statistical information pertaining to the over- 
all population of outflows. In Tables|2]and[3]we list the H2 flows 
catalogued by Sta02 and the newly identified H2 flows from 
this paper, respectively: in both tables we list the likely out- 
flow source (if known), the dense cores that coincide with the 
source, and any HH objects that are associated with the flow. 
In Sect.4. we discuss the region as a whole, drawing statistical 
information from the sample of flows and associated sources. 
In Appendix A we briefly discuss the full sample of molecu- 
lar H2 outflows, presenting continuum-subtracted H2 2. 122 fim 
images of the newly-observed flows. In Appendix B we dis- 
cuss our proper motion measurements for outflows observed 
here and in the earlier work of Sta02. 



2. Observations 

2.1. Near-IR imaging 

Broad-band J, H, K and narrow-band H2 2.122pm images of 
a 1.5° x 1.5° field centred on M 42/M 43 were obtained during 
instrument commissioning on 24 November 2004 (filter char- 
acteristics are described in Hewett et al. 2006). Broad-band K 
and narrow-band H2 images of a ~8 square degree field were 
later secured during Director's Discretionary Time (DDT), on 
15-16 December 2005. K and H2 images of a single 0.75 square 
degree segment centred on the Orion nebula were also obtained 
during Service observing time on 12 February 2007. The com- 
missioning data have only been used to construct the colour 
images in Fig. |2 the later observations, which were secured 
under slightly better observing conditions, and when the in- 
strument was better characterised, have been used for all other 
figures and analysis. The service data have been used to fill in 
a gap in the DDT data, where bright stars and saturation affects 
produced poor results in the processing. 

The near-IR wide-field camera WFCAM (Casal fet al.l 
2007) at UKIRT was used on each occasion. WFCAM houses 
four Rockwell Hawaii-II (HgCdTe 2048x2048) arrays spaced 
by 94% in the focal plane. The pixel scale measures 0.40". 
To observe a contiguous square field on the sky covering 0.75 
square degrees - a WFCAM "tile" - observations at four po- 
sitions are required. At each position, to correct for bad pix- 
els and array artifacts, a five-point jitter pattern was executed 
(with offsets of 3.2" or 6.4"); to fully sample the seeing, at 
each jitter position a 2x2 micro-stepped pattern was also used, 
with the array shifted by odd-integer multiples of half a pixel. 
20 frames were thus obtained at each of the four positions in 
the tile. Eleven tiles in total were observed covering over 8 
square degrees, through both broad-band K and narrow-band 
H2 1-0S(1) filters. Exposure times of 5 sec x 2-coadds and 
20 sec x 2-coadds were used with the K and H2 filters, respec- 
tively. With the jitter pattern and micro-stepping the total on- 
source/per-pixel integration time in K was therefore 200 sec; in 
H2 the total exposure time was 800 sec. 

The DDT data were reduced by the Cambridge 
Astronomical Survey Unit (CASU), which is responsible 
for data processing prior to archiving and distribution by 
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Fig. 1. Our near-IR K-band mosaic of Orion A. In the left-hand panel, the positions of candidate protostars identified from Spitzer 
photometry are marked with circles; the areas outlined with full lines have been mapped at 1200 fj.m (Stanke et al., in prep.); the 
areas marked with dashed lines have been mapped at 850/im ( |Nutter & Ward-Th ompson 2007). In the right-hand panel the 
regions shown in detail in Figs. [3] to [T3] are indicated with boxes. 



the Wide Field Astronomy Unit (WFAU). However, the 
processing of frames that included the bright Orion Nebula 
and Trapezium cluster gave images with severe background 
structure, probably related to the adopted method for sky- 
subtraction. Our service observations of this central region 
were therefore reduced using the ORAC-DR pipeline at the 



telescope (Cavanagh et al. 2003), which performs flat-fielding 
using twilight flats, but does not subtract a sky frame. The 



CASU reduction steps are described in detail by Dye et al. 
(2006). Residual sky/background structure was later removed 
by fitting a coarse surface to each image (as described by 
IDavis et al.l 2007). For both the CASU and Service data, 
astrometric and photometric calibrations were achieved using 
2MASS (|Dye et al.| 2006; IHewett et all 2006); the calibrated 



images were subsequently used to construct the large-scale 
mosaics presented in this paper. 

In our K-band images we reach a limiting magnitude of 
~18.3 in the less-nebulous regions; our sensitivity to point 
sources in the very bright Orion nebula regions is 2-3 magni- 
tudes worse. In the H2 images of L 1641, outflow features with 
a surface brightness of ~ 7xl0~ 19 Wm-2 arcsec~ 2 are detected 



at 3-5cr above the surrounding background. The H2 sensitivity 
in the ONC region is again lower because of the variable, dif- 
fuse nebulosity. 

An overview of the region observed with WFCAM is given 
in Fig.[T] colour images of the spectacular OMC 1/2/3 regions 
are presented in Fig. |2j large-scale H2 1-0S(1) mosaics of re- 
gions of interest throughout Orion A are presented in Figs. [3|- 
[T3l Continuum-subtracted images of the newly-identified H2 
flows are available in Appendix A. In all of these figures axes 
are labelled in 2000.0 coordinates. 



2.2. Proper motion measurements 

Using the newly-acquired WFCAM images and the H2 mo- 
saics of Sta02, which were obtained between December 1996 
and May 1998 with the 3.6-m telescope at the Calar Alto 
Observatory, Spain, we have measured tangential velocities for 
emission-line features in a number of flows. We focus on re- 
gions where multiple flows are observed, particularly where 
the association between outflows and embedded protostars is 
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Fig. 2. Left: colour image of OMC 1/2/3 composed of broad-band J, K and narrow-band Jrb 2.122/zm observations. The data 
have been stretched logarithmically. Right: colour J, H and narrow-band H2 2.122/im composite of the region between Haro 
5a/6a and M 43 (OMC 2/3). The data have again been stretched logarithmically. H 2 flows are labelled with "SMZ" or "DFS" 
numbers; outflow sources are labelled with "IRS" numbers. The dashed box in the left-hand image marks the area covered by the 
right-hand image. 



ambiguous. A complete list of proper motions (PMs) is given 
in Appendix B, where the technique is described in detail. 



2.3. Mid-IR observations 

The Spitzer Space Telescope observations discussed in this 
paper were obtained with the IRAC and MIPS cameras 
dFazio et al.l 2004; lRieke et al.1 2004). A full description of the 
data analysis and young stellar object identification will be pre- 
sented in Megeath et al. (in prep); we provide here a brief sum- 
mary. 

Images in IRAC bands 1, 2, 3 and 4 (at 3.6, 4.5, 5.8 and 
8.0 /zm) and MIPS band 1 (at 24 //m) were used to compile a 
list of more than 300 candidate protostars between declinations 
-4.6° and -9.0°. The protostars were selected from their [3.6]- 



[4.5] and [4.5]-[24] colours fMegeath et al.| 2009), although 
note that the majority of the sources would also be classified 
as flat spectrum, Class and Class I protostars on the basis of 
their Spectral Energy Distributions (SEDs). To minimize con- 
tamination from extragalactic sources we rejected all sources 
with protostellar colours but 24 /mi magnitudes fainter than 7. 
About 20% of the protostars were not detected at 24 /mi and 
were therefore identified through their [4.5]-[5.8] colour using 
the approach of Guterm uth et al.l (2008): most of these sources 



are found in the Orion nebula where the 24 /mi data were satu- 
rated. 

There is a small amount of contamination from extragalac- 
tic sources, particularly AGN; we estimate this contamination 
to be 0.7 sources per sq. degree. In addition, there may also be 
some contamination from embedded Class II sources and po- 
tentially from edge-on young stars with flared disks. Estimates 
of the degree of contamination await the analysis of recent 
Spitzer 5-40 /zm spectroscopy of the sample. In addition to the 
protostars, a further ~ 2000 disk excess sources were identified 
in the Spitzer data. 

The sample of protostars is not complete: we expect that 
many protostars are hidden in the bright nebulosity toward 
the Orion nebula, and that several very bright prototstars were 
omitted because of strongly saturation in the 24 //m data. 

To characterize the SED of each protostar and disk ex- 
cess source, a spectral index was obtained by fitting a single 
power-law to the 3.6, 4.5, 5.8, 8 and 24 /zm photometry. For 
sources without 24 /zm detections, the spectral index was calcu- 
lated from only the IRAC photometry. The resulting 3.6-24 /zm 
spectral-index, a = dlog{AF(A)) / dlogA, ranges from -0.5 to 3.0 
for the protostars. Although the spectral index is useful as an 
indicator of evolution, the index does depend on the density, 
the rotation rate, and the inclination of the protostellar enve- 



lope (Wh itney et al.| 2003, 2004). Consequently, one should not 
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Fig. 3. Narrow-band H2 2. 122 fim images of OMC 2/3, the region directly north of the ONC in Fig.[T](inset figures are double the 
scale and displayed with a logarithmic stretch). Here and elsewhere in this paper Spitzer protostars and disk sources are marked 
with circles and triangles, respectively. MAMBO 1200 /im contours are over-plotted; levels measure 5cr, 7 .5cr, lOcr, 20<r, 40<r, 
80cr, etc., (lcr ~ 20 mJy/beam in central regions of the map where there is no emission). H2 flows are labelled and marked with 
dashed arrows (see also Sta02). 



expect a fool-proof relationship between a and the age of the 
protostar. A histogram showing the distribution of a for the pro- 
tostars in Orion peaks at around and gradually declines with 
increasing a. In contrast, the disk excess sources have spectral 
indices in the range -2.5 < a < 0.0; a histogram of a values 
for these objects peaks at around -0.8, consistent with a large 
population of Class II objects. We note that the spectral indices 
of the disk excess sources can be affected by extinction. Hence, 
the distribution of their a values overlaps the distribution of a 
for the protostars (note that there is not a discontinuity in the 
distribution of a between the disk excess and protostars). This 



may reflect a continuous transition between protostars and disk 
excess sources, as the infalling envelope is dissipated and the 
inner star/disk system is revealed. 



2.4. (Sub)millimetre observations 

In many of the subsequent figures we also overlay contours 
of dust continuum emission at 1200 /mi. These data were ob- 
tained with the IRAM 30-m telescope at Pico Veleta, Spain, 
with the 37 and 117-pixel Max Planck Millimeter Bolometer 



(MAMBO) arrays (Kreysaetal. 1999). The data were ac 
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Fig. 4. Narrow-band H2 2.122/iin images of the region ~20' south of the ONC dJohnstone & Bally (2006) refer to this region as 



OMC 4). Inset figures are again double the scale and displayed with a logarithmic stretch. The full line on the left marks the edge 
of the MAMBO map. DFS 105, DFS 1 10 and DFS 1 14 lie to the west of the region shown above. Contours and symbols are as 
in Fig. [3] 
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Fig. 5. Narrow-band H 2 2.122 fan images of the region around NGC 1980 - OMC 5 in |Johnstone & Bally| (2006) - about 45' 
south of the ONC (see Fig. [T). Inset figures are again double the scale and displayed with a logarithmic stretch. Contours and 
symbols are as in Fig. [3] 



quired over several observing runs between 1999 and 2002, 
and reduced using the in-house reduction package "MOPSI". 
Standard reduction steps were employed; baseline-fitting, de- 
spiking, correction for atmospheric extinction and flux calibra- 
tion using planetary observations. At IRAM the MAMBO Half 
Power Beam Width (HPBW) measures 11". 



A preliminary catalogue of ~500 MAMBO cores was ob- 
tained from visual inspection of the map and 2-dimensional 
Gaussian fitting of features. From this list we identify cores that 
coincide with an H2 outflow source and list the core parameters 
in Tables |2] and |3] The MAMBO dataset will be analysed in de- 
tail in a future paper (Stanke et al., in prep.). 
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3. Results 

3.1. Overview 

In Fig. Q]we present a mosaic of the complete K-band dataset. 
We indicate on this figure the extent of the 850 fim and 1200 fim 
maps, and mark the positions of ~300 protostars. 

Colour images of the bright, nebulous regions around the 
ONC/M 42, the M 43 HII region, and the complex star form- 



ing environment around Haro 5a/6a are shown in Fig. [2] These 
allow the reader to better distinguish H2-emitting jets and out- 
flows from filaments of bright nebulosity. Near-IR observa- 
tions of the region around BN/KL, including the spectacular 
Orion Bullets, the Orion Bar, BN and the Trapezium cluster, 
have been discussed in many other papers (e.g. lUsuda et al.l 
1996; |Schild, Miller & Tennyson] 1997; ILee & Burtonl 2000; 



Bally et al. 2000; Carpenter et al.l 2001; Islesnick et al.l 2004; 
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Fig. 7. Narrow -band H2 2.122yum images of the region around HH 33/40 and HH 34 (roughly 0.5° west of L 1641 -N - see Fig.[TJ. 
The zig-zag line through the centre of the image marks the western border of the MAMBO map. Both insets are again displayed 
at twice the scale and with a logarithmic stretch. Contours and symbols are as in Fig. [3] 



Kas sis et al.l 2006; |Simpson et al.| 2006; ITamura et aU 2006). In 
this paper we therefore focus on the area to the north of the 
ONC, the OMC 2/3 region, and the extensive L 1641 low-mass 
star forming region to the south. 

In Figs.l3lfT3lwe present segments of the full H2 mosaic, in- 
dicating with dashed lines the locations of H2 jets and outflows 
(the data in these figures have been binned over 4x4 pixels, to a 
pixel scale of 0.8", to enhance the fainter features). Many of the 
outflows imaged here have already been catalogued by Sta02 



(H2 images of OMC 2/3 are also presented by lYu et aLl 1997 
and lYu et al.l 2000). With the newly-identified H2 flows, we 
therefore continue the numbering scheme of Sta02, although 
new flows (often found beyond the bounds of the Sta02 data) 
are numbered above 100. The Sta02 objects and the new H2 
outflows are prefixed with SMZ and DFS, respectively. Note 
that entire H2 flows, rather than individual H2 features or small 
groups of knots, are given a single number. For example, HH 1 
and HH 2, which are relatively discrete objects, are here (and in 
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Fig. 8. Narrow-band H 2 2.122;um images of the region around HH 1/2 in NGC 1999 (roughly 0.5° south of L 1641-N - see 
Fig-ID- The inset of HH 1 is double-scale and displayed with a logarithmic stretch. Contours and symbols are as in Fig. [3] 



Sta02) referred to as SMZ 64. All H 2 flows are listed in Tables|2] 
and[3j continuum-subtracted H 2 images of the new DFS flows 
are given in Appendix A, along with a brief description of each 
outflow. 



In Figs. [3lfT3l we also mark the positions of the protostars 
and (more numerous) disk excess sources identified from the 
Spitzer data; candidate outflow sources, when identified, are 
numbered using the same value as is given to the H 2 outflow 
itself (so for example DFS 101 is driven by IRS 101). Outflow 
sources are listed in Tables [2] and [3] together with any coinci- 
dent molecular cores seen at 850 [im or 1200 fj.m. 850 /vm cores 
are taken from |Nutter & Ward-Th ompson (2007); the table of 
1200 fj.m cores has yet to be published (Stanke et al., in prep.), 
so here we again use the outflow number, together with MMS 
(for MilliMetre Source), to identify each core. 



The association of jets with young stars and/or 
(sub)millimetre continuum peaks is largely based on mor- 
phology and/or alignment. However, in some of the more 
complex regions we use the PM measurements in Appendix B 
to help isolate the location of the outflow source; in most cases 
the direction of propagation of the shock features are as one 
might expect, based on their morphology and the location of 
candidate outflow sources. 

3.2. Outflow activity in Orion A 

OMC 2/3 (Figs.|2]and[3l), the region to the north of the ONC, is 
described in considerable detail by lYu et al.l (2000). HH ob- 



jects are identified by Reipurth, Bally & Devine (1997) and 



Ball y & De vine I (2001). High -resolution CO outflow maps are 
presented by [Williams et al. (2003). The distribution of dense 
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Fig.9. Narrow-band H2 2. 122 /mi images of the region ~0.5° east of HH 1/2 in NGC 1999 (see Fig . [TJ . Both insets are displayed 
at double-scale and with a logarithmic stretch. Note that the MAMBO 1200//m (and SCUBA 850 /im) maps do not cover this 
region. Symbols are as in Fig. [3] 



molecular gas throughout OMC 1/2/3 has been mapped at 
450yum and 850/rni by |Johnstone & Bally| (1999); IChini eTahl 
(1997) present 1300z/m observations of OMC 2/3. A compre- 



hensive review is also given by Peterson & Megeath (2008). 

The vast majority of protostars in Fig.[3]are located within 
11" of a MAMBO core and/or 14" of a SCUBA core (these 
radii being equivalent to the HPBW of these (sub)millimetre 
data). The association of bright H2 jets with Spitzer protostars 
and cold molecular cores attests to the youth of the outflows in 
this region. Most of the protostars are located along the chain of 
cores that runs through the centre of OMC 2/3 (the disk excess 
sources are more widely distributed); similarly, most of the H 2 
outflows seem to be driven by sources close to this central axis. 
The most distinctive H 2 flows, SMZ 3, SMZ 5, SMZ 6 and 
SMZ 9 are also orientated perpendicular to the north-south axis 
of cores. 

From Sta02 there are 22 H2 flows spread across the region 
shown in Fig. [3] (we regard SMZ 11 and SMZ 13 as a single 
flow; likewise SMZ 14/16 and SMZ 21/22). We identify four 



new H 2 flows north of the ONC, DFS 101-104. Candidate out- 
flow sources are found for 24 of the 26 H 2 flows. Only two 
have negative spectral indices; the rest are either "flat spec- 
trum" or reddened protostars. Two of the 25 Spitzer outflow 
sources were not observed with MAMBO; of the remaining 
23, 18 are associated with 1200 /zm (MMS) cores. 

In Figs. |4]and[5]we show large-scale H 2 images of two re- 
gions south of the ONC, adjacent to NGC 1980 in the overview 
plot in Fig. Q] Johnstone & Bally (2006) refer to these regions 
as OMC 4 and OMC 5, respectively. In and around these two 
areas Sta02 identify at least 18 possible H 2 flows (SMZ 29- 
SMZ 46; between declination -5°20' and -6° 15'); we add a 
further ten H 2 flows to this tally (DFS 105-114; note that 
DFS 105, DFS 110 and DFS 114 are beyond the edges of 
Fig. |4j. Moreover, PMs have been measured for many of the 
flows in this area. 

We identify candidate outflow sources for 25 of the 28 out- 
flows in the OMC 4 and 5 regions, although unlike OMC 2/3, 
many of these sources are not associated with dust cores. (Note 
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Fig. 10. Narrow-band H2 2.122/im images of L 1641-C (labelled in Fig. [TJ. Both insets are displayed with a logarithmic scale 
and at twice the scale of the main figure. The northeast and northwest edges of the MAMBO map are marked with full lines. 
Contours and symbols are as in Fig. [3] 



that Johnstone & Bally (2006) also compare H2 images with 
dust continuum maps in this region, using the images of Sta02 
and their SCUBA observations, respectively.) All of the SMZ 
jets were observed with MAMBO, as were seven of the 10 DFS 
flows. However, only 14 of the 25 observed H2 outflow sources 
are associated with 1200 /im cores. There are clearly consider- 
able differences between the star forming regions north of the 
ONC (OMC 2/3), and those directly south (OMC 4/5). 

As we move further south (through L 1641-N towards 
L 1641-C, see Fig. [TJ fewer protostars are identified in the 
Spitzer analysis. Even so, many sources are clearly driving 
H2 flows: in many regions (e.g. around L 1641-N, HH 34 and 
HH 1/2 in Figs. [6] Hand [8} we note a clear association between 
Spitzer protostars, (sub)millimetre cores and H2 outflows. 

The infrared jets and outflows in the very busy re- 
gion around L 1641-N (Figj6j are described in detail by 



iGalfalk & Olofssonl (2007) and IStanke & Wifliamsl (2007). 
The former present Spitzer observations and PMs for knots 
close to L 1641-N itself; the latter show high-resolution CO 
maps that reveal a number of bipolar outflows centred on 
this cluster. The HH objects in the region are discussed 



by |Reipurfh, Devine & Bally| (1998). Sta02 identify seven H 2 
flows in and around L 1641-N (SMZ 48-SMZ 54); we label 
two additional flows, DFS 118 - which is associated with a 
bipolar CO outflow ( IStanke & Williamsl 2007), and DFS 115 
- which is associated with a small cluster of protostars ~3' 
south of L 1641-N, known as Strom 11. Six of these nine H2 
flows have candidate Spitzer protostellar sources; all but one 
of these protostars is associated with 1200/im emission or a 
MAMBO core. In total there are 10 protostars within a 5' ra- 
dius of L 1641-N; it seems likely that most, if not all of these 
sources could be associated with H2 line emission features. 
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Fig. 11. Narrow-band H2 2. 122yum images of the region ~0.3° west-southwest of the L 1641-C cluster that includes the extensive 
HH 38/43/64 bipolar outflow (see Fig. [T). HH 43 and its source are shown inset at double-scale with a log stretch. The northern 
and western edges of the MAMBO map are marked with a full line. Contours and symbols are as in Fig. [3] 



The region to the west of L 1641-N, around HH 34, 
is shown in Fig|7j The source of the well-known HH 34 
jet, HH 34-IRS (IRS 55), is flagged as a protostar in the 
Spitzer data, along with two other sources, one of which 
(IRS 56) drives a second H 2 flow. North of HH 34, the 
well-known HH objects HH 33/40 have in recent years been 
considered to be part of a parsec-scale flow that includes 
HH 34 and HH 8 6/87/88 to the south (IDevine et al.l 1997; 
Eisloffel & Mundt 1997). However, in the Spitzer analysis a 
protostar (IRS 1 16) is found coincident with HH 40. The source 
is point-like in 24 /mi MIPS data so is probably not a mis- 
identified H2 feature^ IRS 1 16 is not associated with a dense 



1 H 2 features produce pure-rotational emission lines in each of the 
four Spitzer IRAC bands, though not in the MIPS 24pm band (e.g. 



ISmith&Rosenl 2005: Velusamy, Langer & Marsh| 2007). Spectral in- 
dices measured from the IRAC bands alone can therefore mimic val- 



SCUBA or MAMBO core. Even so, it remains a viable candi- 
date outflow source for HH 33/40. 



Roughly 25' west of HH 34, HH 83/84 ( Reipurth 1989) 
is also observed in H2 emission. HH 83-IRS, the source of 
this well-known jet, is outside the bounds of the SCUBA 
and MAMBO observations, and the Spitzer IRAC 3.6 yum and 
5.8 jt/m and MIPS 24yum maps, although a bright source is de- 
tected in the remaining IRAC bands with an infrared excess. 
HH 83-IRS's [4.5]-[8] color is more consistent with a disk ex- 
cess source than a protostar. 

The outflows in the vicinity of HH 1/2, shown here in Fig. [8] 
are well documented (see e.g. Davis, Eisloffel & Ray 1994; 

Reipurth et al. 



Hester, Stapelfeldt & Scowen 



Radial velocities and PMs 



1998; |Reipurth et al] 2000) 
based on optical observations - 



ues expected for protostars (Davis et al. 2008), though a MIPS detec- 
tion would then require bright, forbidden line emission from [Fell]. 
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Fig. 12. Narrow -band H2 2.122yum images of the region around Re 50 and Haro 4-255, midway between L 1641 -C and L 1641 -S 
in Fig. 03 The figures inset show the regions around Haro 4-255 and the likely source of SMZ 76 (both displayed with a log 
stretch/double size). The southern edges of the MAMBO map are marked with a full line. Contours and symbols are as in Fig. [3] 
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Fig. 13. Narrow-band H2 2.122yum images of L 1641-S, at the bottom of the overview mosaics in Fig. Q] No MAMBO (or 
SCUBA) observations are available for this region. The figures shown inset are displayed at twice the scale and with a log 
stretch. Symbols are as in Fig. [3] 



have been reported in the literature for the central bright HH 
flows dEisloffel, Mundt & Bohm] 1994; [Bally et al.| 2002), so 
we have only measured PMs for the bright H2 features 4'- 
5' northeast of HH 1/2. We find that these knots are moving 
southward, which confirms their proposed association with the 
large-scale SMZ 49 outflow. Around HH 1/2 itself, only four 
Spitzer protostars are identified (SMZ 62/HH 3 is probably a 
mis-identified H2 knot). Three of the four protostars are associ- 
ated with 1200 /im cores; all four sources likely drive H2 flows. 
The well-known VLA source of the HH 1/2 outflow is unre- 
solved from neighbouring sources and therefore is not listed as 
an "IRS" source in Table [2] 

By extending the maps to the southeast and southwest of 
the HH 1/2 region we discover six more H2 flows, DFS 120 
and DFS 122-DFS 126, four of which are shown in Fig. [9] 
(DFS 122 is associated with HH 89; |Reipurth| 1985). Three 
of the six have candidate protostellar sources (one was not ob- 
served with Spitzer). With a length of 1.1 pc, DFS 123 qualifies 
as a "parsec-scale" flow. 

Roughly 30' to the southeast of HH 1/2, the Spitzer data 
recover the cluster of young stars collectively known as L 1641- 
C (Fig. [Toll. However, like Sta02 we identify only one clear-cut 



flow in the region, SMZ 66. The cluster is associated with only 
faint 1200 fim emission and so L 1641-C probably harbours 
a fairly evolved population. The absence of ambient material 
explains the lack of H2 flows. 

FigCEUshows the parsec-scale flow HH 38/43/64 (SMZ 67). 
The source of this outflow, and the neighbouring H2 jet 
SMZ 68, are both retrieved from the Spitzer data. These pro- 
tostars are also associated with molecular cores. If we include 
the DFS 126 bow-shock in the HH 38/43/64 outflow, then the 
total length of this system is 25' (3.2 pc). 

The region midway between L 1641-C and L 1641-S 
(FigfT2l see also Fig. [JJfor reference) includes the well-known 



T Tauri star Haro 4-255 ( |Aspin & Reipurth 2000), its infrared 
neighbour Haro 4-255 FIR, and the spectacular arc of H2 emis- 
sion, SMZ 76, first noted bv lStankel (2000). SMZ 74 (not shown 
here), SMZ 75 and SMZ 76 all have morphologies and PMs 
that imply motions away from the chain of young stars near the 
centre of FigQ~2l although precise associations are difficult to 
make. [Stankel (2000) identified a millimetre peak southwest of 
SMZ 76 (which he labeled L 1641-S3 MMS1) as the location 
of the likely source of this object; we note that DFS 132 may be 
the counter-lobe of this remarkable outflow. The protostar/dust 
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Fig. 14. Histograms showing (a) the distribution of jet lengths 
(in 1 ' bins), (b) jet opening angles (5° bins), and (c) jet position 
angles (10° bins), for all H2 jets in Orion A (open columns), for 
H 2 flows north of the ONC (declination > -5°30'; grey filled 
columns) and for H 2 flows south of the ONC (< -5°30'; black 
filled columns). 



core IRS 131/MMS 131 may also drive an H 2 flow (DFS 131), 
though there are perhaps a dozen Spitzer YSOs in this region 
with no obvious H 2 jet. Most have no associated dust core and 
are therefore, like the young stars in L 1641-C, probably too 
evolved. 

Lastly, the region at the southern end of L 1641 is abun- 
dant with Spitzer-identified protostars and newly discovered H 2 
flows. In Fig[13]we label six possible outflows, four of which 
have candidate sources (note that the MAMBO observations do 
not extend this far south in L 1641). DFS 138 is the most spec- 
tacular H 2 flow, comprising a bright, collimated, bipolar jet and 
sweeping bow shocks. The flow extends over 13.4' (1.8 pc) and 
is probably driven by the Spitzer protostar labelled IRS 138. 



Fig. 15. A map showing the distribution of short (L <1'; open 
squares), intermediate (l'< L < 2'; filled grey squares), and 
long (L >2'; filled black circles) H 2 jets. 



3.3. Outflow statistics 

In Tables[2]and[3]we list 1 16 jets, of which 43 are newly iden- 
tified. Of the 111 flows that were within the bounds of the 
Spitzer observations, 72 (62%) have catalogued Spitzer proto- 
stellar sources: a further 12 flows are tentatively associated with 
Spitzer sources (IRS source numbers marked with a question 
mark in Tables [2] and |3j. Two well-defined flows, HH 1/2 and 
V380 Ori-NE, have established sources that were not flagged as 
young stars in the analysis of the Spitzer data (the former be- 
cause it is too faint; the latter because it is too close to a neigh- 
bouring source); a third object, HH 83, was only observed in 
two of the four Spitzer IRAC bands. Closer examination of the 
mid-IR data shows that HH 1/2 and V380 Ori-NE are probably 
driven by protostars, and HH 83 by a disk excess source. 

In total, 97 H 2 outflows (including all of the SMZ jets) were 
observed with MAMBO at 1200 /zm. 53 (55%) are associated 
with MAMBO (MMS) cores and young Spitzer sources; a fur- 
ther four outflow sources are associated with extended 1200 jt/m 
emission or are unresolved from brighter, nearby cores. 

Of the 116 jets, 68 (59%) exceed an arcminute in length 
(at a distance of 450 pc, 1' is equivalent to 0.13 pc); 10 are 
parsec-scale jets, exceeding 7.6' (SMZ sources 6, 49, 55, 67, 
76 and the newly-discovered DFS flows 117, 123, 132, 133 and 
138). Throughout Orion A we measure a range of jet lengths (as 
traced in H 2 ); a histogram showing the distribution of lengths, 
uncorrected for inclination angle, shows that the number of 
flows with a given length decreases more-or-less exponentially 
(Fig. fT4h). The same trend is seen if we just consider flows 
north of the ONC or flows south of the ONC: the map in Fig. [TBI 
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Fig. 16. Source spectral index, a, integrated 1200 fim core flux, /1200, and 1200 fim flux density, F1200, plotted against H2 jet 
length and jet opening angle for all jets in Orion A with identified sources (Tables[2]and[3]). 



shows that the distribution of short (L < 1'), medium-length, 
and long (L > 2') jets is more-or-less even across Orion A. 

The jet lengths in Tables [2] and [3] are of course lower lim- 
its, since they are uncorrected for inclination to the line of sight. 
Also, as a jet or wind exits its molecular cloud, it will no longer 
entrain and shock-excite molecules and may not be visible in 
H2 emission. The histograms in Fig. [T4h are therefore biased 
towards shorter flow lengths. Even so, if all of the jets were 
of equal length and randomly inclined with respect to the line 
of sight, inclination alone would produce a histogram that in- 
creases towards longer jet lengths (peaking at an angle of n/2, 
neglecting any selection effects due to detectability as a func- 
tion of angle). We can therefore be sure that the "H2 flows" in 
Orion are mostly short; macro-jets are a rare occurrence and, 
statistically at least, H2 jet lengths do not appear to be dictated 
by the changing environment in Orion A. 

We find no real correlation between H2 jet length and indi- 
cators of the outflow source age, such as the spectral index, a 
(defined in Sect. 2.3), the integrated 1200 yum flux of the core 
associated with each outflow source, /1200, or the 1200 yum sur- 
face brightness measured towards each outflow source, F1200 
(Fig. Q2]- note that F1200 is not always the flux density mea- 
sured towards the core peak). Five of the six longest H2 flows 
are associated with IRS sources with large spectral indices 
(a > 1 .5), though at the same time four of these five protostars 
seem to be associated with modest-sized cores (/1200 < 1-5 Jy; 
^1200 < 0.6 Jy) - the fifth protostar was not observed with 
MAMBO. The longest H 2 jet in Orion A, SMZ 49 (L ~ 35'), 
is almost certainly driven by a source in the compact L 1641- 
N cluster, which is associated with highly-reddened protostars 
and dense, massive molecular cores. However, identifying this 
remarkable jet's progenitor and, in particular, resolving the as- 
sociated core, is difficult. 



There is some indication from millimetre observations 
of CO outflows that Class sources produce more colli- 
mated flows than their Class I counterparts dLee et al.l 2002; 
|Arce & Sargent| 2006). This may be due to the destruction of 
the protostellar core by the outflow itself, as it entrains and 
plows away the ambient gas. Since H2 emission traces the in- 
teraction of outflows with the surrounding molecular gas, one 
might expect to see some indication of this "flow broadening" 
in H2 images of outflows from sources with lower values of a, 
or low values of 7i2oo- However, this does not seem to be the 
case: the flow opening angle, defined as a cone that encloses all 
H2 emission line features, with an apex centred on the outflow 
source, appears to be unrelated to a, Inoo or F1200 (Fig. [Toll. 

If we simply compare sources associated with a molecular 
core with sources undetected in 1200/im emission, we again 
find that both samples essentially have the same distributions 
of H2 jet lengths and jet opening angles: for H2 jet sources 
with cores the mean length is 3.9' (standard deviation = 5.6'; 
range = 0.2'-34'); for those without cores the mean length is 
1.7' (std dev = 1.9'; range = 0.1 '-9. 5'). The former is skewed 
somewhat by the fact that four of the seven longest H2 flows 
(L >10') are associated with IRS sources that do coincide with 
1200 /mi cores (the other three long flows were not observed at 
1200 fim). The mean opening angles for H2 jets with and with- 
out cores are 27.7° (std dev = 18°; range = 4°-70°) and 19.9° 
(std dev =15°; range = 3°-50°), respectively. 

The lack of a correlation between a and H2 jet length, L, or 
opening angle, 8, may be a result of the fact that a varies with 



orientation as well as source "age" (e.g. Whitney et al. 2003). 
More evolved sources driving the longest jets are expected to 
have lower spectral indices. However, jets inclined towards the 
observer will appear shorter and will likewise have lower spec- 
tral indices. Poor resolution and source confusion may also ex- 
plain the lack of a correlation between core parameters (/1200 
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and F1200) and L or 6. One should also note that H2 is not al- 
ways an ideal measure of flow parameters. As noted earlier, 
H2 excitation only occurs when the flow impacts dense ambi- 
ent material, and pre-existing H2 features cool and fade very 
rapidly (within a few years). For a flow that exits a molecular 
cloud, H2 emission features will not be observed, so the true jet 
length may be underestimated. Conversely, jet opening angles 
may be over-estimated: the wings of jet-driven bow shocks are 
often much wider than the underlying jet (SMZ 42 is a nice 
example of this in Orion), and changes in flow direction due to 
precession will, over time, increase the apparent opening angle. 

Finally, the precise relationship between a and source age 
has not yet been established. Spectral typing of protostars is dif- 
ficult because of extinction and veiling of photospheric absorp- 
tion lines (IGreene & Ladal 2002: Doppmann et al. 2003; lBeckl 
2007; IAntoniucci et al.1 2008). so age estimates based on the fit- 
ting of theoretical mass tracks and isocrones, especially for the 



youngest sources, can be problematic (see e.g. Doppmann et al. 



2005). Radio continuum observations may be a better way to 
quantify the youth of an outflow source. 

3.4. Flow orientations 

H2 jet position angles on the sky should be well defined, 
provided the outflow source is correctly identified. H2 sur- 
veys can therefore be used to test the standard paradigm that 
clouds collapse along field lines to form elongated, clumpy fil- 
aments from which chains of protostars are born (Mouschovias 
1976; [Larry et al.| 2008), with the associated outflows aligned 
parallel with the local field, and perpendicular to the chains 
of cores (Banerjee & Pudritz 2006). The presence of mag- 
netic fields in molecular clouds has been demonstrated, and 
the orientations of field lines with respect to chains of oblate 
cores have been measured in a number of star forming regions 
(e.g. ICrutcheret all 19 99; IMatthews etaH 2002; IHoude et al.l 
2004; |Vallee & Fiege] 2006; |Larry et al.| 2008). In a few re- 
gions there is also some indication that outflows are prefer- 
entially orientated orthogonal to individual cores or chains of 
cores, and/or parallel with the local magnetic field dStrom et al.l 
1986;|Heyer, Strom & Strom] 1987; |M£ri^d_&puchenel 2004; 
IDavis et al.l 2007; |Ananthapiridika & Whitworth| 2008). 

The Orion A GMC comprises a large number of cores 
extended along a north-south "integral-shaped" filament that 
passes through the Orion nebula and extends southward, broad- 



ening out between L 1641-N and L 1641-S (Johnstone & Bally 
1999, 2006; |Nutter & Ward-Thompson] 2007). It is therefore 



potentially an ideal laboratory for an investigation into the re- 
lationship between cores, filaments and outflows. Throughout 
much of Orion A, one might expect to see flows predominantly 
orientated east-west. 

In Fig. fT4b we have plotted a histogram of the jet posi- 
tion angles (in ten degree bins) for all outflows in our sam- 
ple. The plot suggests a uniform distribution of position angles. 
Moreover, a Kolmogorov-Smirnov (KS) test shows that there is 
a >99% probability that such a distribution is drawn from a ho- 
mogeneously distributed sample. In other words, the position 
angles of the H2 flows are completely random. If we consider 



only the jets with identified outflow sources a similar situation 
is found. In this case the probability that the distribution is ran- 
dom is 99.2%. The same also holds if we separate the H2 jets 
in the northern region (Dec > -6°), where the dense gas is more 
tightly constrained to a north-south filament, from those south 
of the ONC (Dec < -6°). Even for the compact ridge of pre- and 
proto-stellar cores north of the ONC (OMC 2/3), the outflows 
are randomly orientated. Similar plots for H2 jets exceeding an 
arcminute in length (flows that arguably have more reliable po- 
sition angles), or for jets with embedded protostars coincident 
with MAMBO cores, are equally randomly orientated. 



So the H2 jets in Orion A appear to be randomly orientated. 
Notably, IMatthews etaD (2001) and IHoude et al.l (2004) find 
that on large scales the field along the integral-shaped filament 
(north and south of the ONC) is non-uniform. The long axes 
of the elongated cores along the filament also appear to be ran- 
domly orientated dChini et al.l 1997; |Johnstone & Bally| 1999, 
2006;|Nutter & Ward-Thompson|2007; see also e.g. Fig.©. It 



is therefore not too surprisingly that flow position angles should 
also be non-uniform. But are the H2 outflows orientated or- 
thogonal to the cores associated with each outflow source? At 
450 pc, the MAMBO beam of 11" corresponds to 5 x 10 3 AU, 
considerably larger than accretion disk dimensions or the size 
of the solar system. Even so, many of the cores in Tables|2]and 
[3]are elongated; if these "circumstellar envelopes" are magnet- 
ically tied to the accretion disks they harbour (due to modest 
ambipolar diffusion), then one might expect them to be orien- 



tated perpendicular to the flow axes (Banerjee & Pudritz 2006; 



Ananthapindika & Whitworth 2008). 



In our data, 47 H2 jet sources are associated with 1200 /mi 
emission above a surface density F1200 > 75 mJy. Of these, 42 
are associated with cores (rather than "emission" in Tables [2] 
and |3). From this sample we extract sources where the core 
major- to-minor axis ratio exceeds 1.3, and exclude very large 
(diffuse) cores where the Full Width Half Maximum of the 
Gaussian fit is greater than 50" in either axis. We arrive at 
a sample of just 22 H2 flows with cores. Surprisingly, even 
from this finely-tuned sample, the distribution of H2 jet posi- 
tion angles (PAs) with respect to core PAs is completely ran- 
dom. One might expect that the "longer", more clearly-defined 
jets, would be orthogonal to their associated cores. However, in 
a plot of jet length against "jet-PA minus core-PA" (not shown), 
even the longest jets seem to be randomly orientated with re- 
spect to their progenitor cores. 



It seems likely that, in some of the more clustered re- 
gions (OMC 2/3, L1641-N, etc.), the MAMBO observations 
have insufficient resolution to disentangle the molecular en- 
velopes associated with multiple sources. Focusing on a hand- 
picked sample of sources separate from the main clusters of 
protostars, sources with well-defined jets and bright, compact 
cores, is probably necessary. Obviously much higher resolution 
(sub)mm observations would also be beneficial; ALMA will be 
useful in this respect. 
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Fig. 17. Histograms showing the number of disk excess sources 
(open columns), protostars (filled grey columns) and H2 flow 
sources (black columns) with a given value of a (in bins of 
0.2). 



emission. Thus, at first sight it does not appear that all proto- 
stars drive H2 flows. 

However, our near-IR images may not be sensitive enough 
to detect the weakest flows, especially H2 jets from the least lu- 
minous - and least massive - protostars dCaratti O Garatti et al.l 
2008). Also, the selection criteria used for identifying the 
Spitzer protostars may well include a handful of background 
galaxies and AGN (~0.7 per square degree), H2 emission knots, 
and a number of more evolved young stars. Though faint, heav- 
ily reddened Class II sources and disk excess sources viewed 



edge-on can have the same colours as protostars (Whitney et al. 



2003). These classical T Tauri stars rarely produce extensive H2 
flows, because they have either swept away much of their cir- 
cumstellar environment, or because they occupy regions with 
a low ambient molecular gas density (they may of course 
drive HH jets, which our near-IR survey is not sensitive to). 
IGutermuth et al.l (2009) estimate that at most 3.6% of proto- 
stars may be edge-on disk excess sources. Our WFCAM im- 
ages include 1909 disk sources (Megeath et al., in prep.). There 
may therefore be as many of 70 disk sources in the protostar 
catalogue that are too evolved to drive H2 outflows. 



4. Discussion 

4.1. Do all protostars drive H 2 outflows? 

In Fig. [17] we plot histograms showing the distribution of spec- 
tral indices for the Spitzer disk excess sources, the Spitzer pro- 
tostars, and the outflow sources. As was found to be the case in 
Perseus (Davi s~etaLl 2008) and the more distant, high mass star 
forming region, W 75/DR 21 (IDavis et al.l 2007: iKumar et al.l 
2007) the majority of H2 flows are driven by protostars with 
positive spectral indices. The mean value of a for the disk ex- 
cess sources, most of which will be T Tauri stars, is -0.85 (std 
dev = 0.61); for the protostars and the outflow sources, the 
mean values are 0.80 (std dev = 0.72) and 0.86 (std dev = 0.93); 
in terms of spectral index, the outflow sources in Orion A ap- 
pear to be indistinguishable from the protostellar sample. 

For the outflow sources in Perseus and DR 21/W 75 the 
mean values of a are 1.4 and 1.9, respectively, somewhat 
higher than in Orion A. The protostellar population in Orion A 
may therefore be slightly more evolved than in Perseus and 
DR 21/W 75. Note, however, that the DR 21/W 75 spectral in- 
dices were derived from only IRAC photometry, which may re- 
sult in a higher mean value (the Perseus indices were measured 
using 2MASS JHK photometry, as well as Spitzer and MIPS 
data). Also, the larger UV radiation field from the ONC and 
from B- and A-type stars spread throughout Orion A may ac- 
celerate the erosion of protostellar envelopes, thereby reducing 
measured values of a, although a is largely defined by emis- 
sion from the inner disk rather than the outer, cold envelope, so 
this affect may be modest. 

The spectral index is clearly a useful tool when search- 
ing for H2 outflow source candidates. Yet of the 290 Spitzer- 
identified protostars observed in Orion A with WFCAM (ex- 
cluding the dozen sources in the Orion Nebula), only ~84, or 
less than one third, seem to be driving outflows that excite H2 



The fact that the distributions of a values for the out- 
flow sources and protostars in Fig. [T7] are comparable sug- 
gests that protostars and H2 outflow sources are essentially the 
same objects; as noted earlier, the spread in a for both sam- 
ples may be due to inclination angle as much as to source age 
( [Whitney et al] 2003; lAUen et al.l 2004; |Megeath et aL] 2004). 
Moreover, there is empirical evidence which suggests that the 
H2 outflow phase may be almost as long-lived as the proto- 
stellar phase: the dynamical age of an H2 flow is of the or- 
der of 10 4 -10 5 years (based on a canonical jet length of 1 pc 
and molecular outflow speeds of 10-100 km s _1 ; e.g. lArce et al.l 
2007), while statistical studies indicate that the low mass pro- 
tostellar phase may last at least as long (Froebr ich et al.1 2006: 
IHatchell et all 2007: Enoch et al.1 2008). 



Towards the end of the protostellar phase the accretion rate 
- and therefore the mass loss rate - will evolve quite dramat- 
ically, as the source approaches the Class II T Tauri phase. 
Statistically at least, these sources will probably be associated 
with mostly atomic, low power jets, which may then persist for 
timescales of the order of 10 5 - 10 6 years (Reipurth & Bally 



2001). Such jets may be undetectable without very deep H2 
imaging, though their sources will in any case usually be iden- 
tified as "disk excess sources" in the Spitzer analysis used here. 

In conclusion, it seems likely that for much of their evo- 
lution all protostars will produce outflows bright in Hi line 
emission. In a sample such as ours, the protostars not associ- 
ated with H2 flows may be misclassified disk excess sources, 
AGN or even H2 knots. Some protostars may be marginally 
too evolved or too low mass to produce detectable H2 jets. 
Orientation effects and the local, inhomogeneous environment 
will also play a role in jet detectability. 
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4.2. Association of H 2 outflow sources with molecular 
cores 



MAMBO 1200 micron Image 



In a recent study, [J0rgensen et al.| (2007) used a combination 
of Spitzer IRAC and MIPS photometry and 850/im SCUBA 
observations to identify the youngest protostars in Perseus. 
They selected 49 MIPS sources that have red [3.6]-[4.5] and/or 
[8.0]-[24] colours, and are within 15" of a SCUBA peak. 
IDavis et all (2008) later found that most (63% of the reddened 
MIPS sources they imaged) drive H2 flows. In Orion A it seems 
that a similarly high percentage of "protostars with cores" drive 
molecular (H2) outflows; note that most of the outflows with 
identified progenitors in Tables |2]and[3] are associated with dis- 
crete 1200 //m cores. 

However, clearly not all molecular cores are associated 
with H2 flows. We identify in excess of 500 peaks in the 
1200 pm MAMBO data (Stanke et al., in prep.). Similarly, 



Nutter & Ward-Thompson (2007) extract almost 400 cores 



from their more limited 850 pm SCUBA maps of Orion A. 
Only about l/5th of these cores seem to be associated with H2 



outflows. Johnstone & Bally (2006) obtain a similar result in 
their comparison of SCUBA data with the H2 images of Sta02; 
they find that 17 of the 70 SCUBA clumps they identify are 
clearly associated with H2 emission. 

Many of the cores in Orion A will be prestellar - based on 
the number of cores without H2 jets, roughly 80% - which sug- 
gests that the H2 outflow (and therefore protostellar) phase is 
considerably shorter than the prestellar phase. This statistic as- 
sumes that each core produces just one protostar (and H2 flow), 
which is almost certainly incorrect; small groups of protostars 
are often found in each core (see for example the L 1641-N 
core in Fig.[6]l. The prestellar phase could therefore be at least 
an order of magnitude longer than the protostellar/H2 outflow 
phase. 

What fraction of cores containing protostars drive H2 
flows? In Fig. [18] we plot the entire 1200 pm map of Orion A. 
We over-plot the positions of the Spitzer candidate protostars, 
and mark the locations of the H2 outflow sources. (Note that the 
outflow and protostar catalogues are incomplete in the Orion 
Nebula region; saturation effects, source confusion, and the 
bright, diffuse nebulosity that pervades the ONC will limit both 
samples.) Of the 190 protostars in Orion A that were observed 
with MAMBO, excluding the dozen in the Orion Nebula re- 
gion, only 70 coincide with 1200yimi cores (to within an 11" 
radius), although most do coincide with extended and/or diffuse 
emission. Of these 70 protostars with cores, 43 (61%) appear 
to be driving molecular H2 outflows. In some of the less tightly 
clustered and therefore less complex regions in L 1641, such as 
the NGC 1980 region in Fig. [5] the number of protostars associ- 
ated with cores that also drive H2 jets approaches 100%. Within 
the limitations of our millimetre and near-IR observations (dis- 
cussed in the previous section), it seems likely that most, if not 
all protostars within molecular cores drive Hi outflows. It is 
also clear that the combination of Spitzer mid-IR photometry 
and either SCUBA 850/ari or MAMBO 1200 pm far-IR pho- 
tometry is ideal for identifying molecular outflow sources and, 
in all probability, the youngest protostars. 
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Fig. 18. A logarithmically-scaled image of 1200 pm emission 
with, superimposed, the locations of all Spitzer-identified pro- 
tostars (open circles). H2 outflow sources are marked with 
crosses (note that a handful of H2 jet sources are disk excess 
sources). 



Lastly, we mention that, of the 30 or so protostars in Fig. [18] 
that are located to the east or west of the edges of the MAMBO 
map, i.e. beyond the bounds of the high-density molecular 
ridge that constitutes OMC 2/3 and L 1641, only three proto- 
stars are driving H2 flows. A much larger fraction of protostars 
within the bounds of the MAMBO map, where ambient gas 
densities are much higher, are driving molecular outflows. As 
was noted in the previous section, some of the protostars with- 
out H2 jets may be marginally more evolved; Fig. [T8l suggests 
that these "evolved protostars", sources close to the Class II 
T Tauri phase (which may of course drive HH jets), are more 
widely distributed. 

4.3. Clustering and the distribution of protostars and 
outflow sources. 

Qualitatively speaking, if one examines Fig. [18] closely, it ap- 
pears that clusters of protostars are not always associated 
with clusters of H2 outflows. In L 1641-C, for example (RA 
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Table 1. Mean and median Nearest Neighbour (NN) radii for protostars (all sources, sources with H2 outflows, and sources 
without) in Orion A. The standard deviation (std) is given in brackets. North and South refer to protostars located north or south 
of declination -6°, respectively. 



Region 


5th NN 


6th NN 


7th NN 


10th NN 


20th NN 




mean(std)median 


mean(std)median 


mean(std)median 


mean(std)median 


mean(std)median 



All Protostars 
Protostars with flow 
Protostars no flow 



0.087°(0.09°)0.060° 
0.077° (0.06° )0.064° 
0.090° (0.10°)0.055° 



0.097°(0.10°)0.070° 
0.089°(0.07°)0.085° 
0.100°(0.11°)0.061° 



0.106°(0.10°)0.079° 
0.098° (0.07° )0.096° 
0.109°(0.11°)0.073° 



0.13°(0.11°)0.11° 
0.13°(0.08°)0.12° 
0.13°(0.12°)0.09° 



0.20°(0.11°)0.18° 
0.19°(0.10°)0.20° 
0.20°(0.12°)0.17° 



All Protostars (North) 
Protostars with flow (North) 
Protostars no flow (North) 



0.077° (0.08° )0.042° 
0.059°(0.06°)0.035° 
0.083°(0.09°)0.047° 



0.086° (0.09° )0.050° 
0.068°(0.06°)0.048° 
0.093°(0.10°)0.051° 



0.095°(0.10°)0.058° 
0.075°(0.07°)0.058° 
0.103°(0.11°)0.059° 



0.12°(0.11°)0.08° 
0.10°(0.08°)0.11° 
0.12°(0.12°)0.08° 



0.16°(0.12°)0.14° 
0.14°(0.10°)0.14° 
0.17°(0.13°)0.13° 



All Protostars (South) 
Protostars with flow (South) 
Protostars no flow (South) 



0.096°(0.10°)0.075° 
0.094°(0.06°)0.093° 
0.097°(0.11°)0.073° 



0.107°(0.10°)0.094° 
0.110°(0.06°)0.120° 
0.106°(0.11°)0.079° 



0.115°(0.10°)0.101° 
0.119°(0.06°)0.137° 
0.114°(0.11°)0.083° 



0.14°(0.11°)0.12° 
0.16°(0.07°)0.18° 
0.14°(0.12°)0.11° 



0.23°(0.11°)0.21° 
0.24°(0.06°)0.23° 
0.23°(0.12°)0.19° 



~ 5*38.7 m , Dec ~ 7°00'; see also Fig.QJJl, although a dozen 
protostars and ~30 disk excess sources are identified, there 
is only one H2 flow. Likewise, between Re 50 and Haro 4- 
255 (RA ~ 5*39.9 m , Dec ~ -7°28'; see also Fig. [H there 
is a chain of 16 protostars and at least as many disk excess 
sources, yet again we only detect three or four molecular H2 
outflows. Neither region is associated with large columns of 
ambient molecular material (L 1641-C is barely detected in the 
MAMBO data); the paucity of detected outflows may there- 
fore be due to the evolved nature of the protostars in each re- 
gion. L 1641-N, on the other hand, is associated with massive 
dense cores, and a large fraction (roughly 50%) of the proto- 
stars power H2 outflows. Likewise, north of the ONC, in the 
OMC 2/3 molecular cloud region (between declinations -4°30' 
and -5° 15') where protostars and cores are very highly clus- 
tered, there are 65 protostars of which 26 appear to be driving 
molecular flows (this may be a lower limit, given the abun- 
dance of H2 emission line features, although protostars are also 
difficult to resolve in this busy, nebulous region). 

We also note that regions with multiple H2 flows often 
seem to be more sparsely populated with protostars. For ex- 
ample, around HH 1/2 (RA ~ 5 ft 36.3 m , Dec ~ -6°46'; see also 
Fig. [8]) there are at least half-a-dozen H2 jets, yet only about the 
same number of protostars. Similarly, east of HH 1/2 around 
DFS 123 (Fig. 13 and south of HH 1/2 around HH 38/43/64 
(Fig. [Toll the ratio of H2 jets to protostars is roughly one to one, 
though the number of protostars is again low. 

However, neither of the above semi-quantitative observa- 
tions seems to be related to clustering. We demonstrate this by 
measuring the radius to the 5th, 6th, 7th, 10th and 20th nearest 
neighbour (NN) for each of the protostars identified in Orion 
A. The distribution of 6th NN radii is displayed in Fig. [19] 
where we distinguish protostars that do drive H2 flows from 
those that do not. We also separate protostars north of Dec -6° 
(the OMC 2/3 region) from those south of Dec -6° (L 1641). 
In Table[JJwe list the mean and median radius for each sample 
of Nth nearest neighbour radii, and again distinguish outflow 
sources from the remaining sample of protostars, and contrast 
the northern region with the south (median values tend to ex- 



clude outliers with large radii located on the periphery of each 
region). Table [TJ and Fig.[l9]both demonstrate that there is no 
significant difference between the spatial distribution of proto- 
stars that drive H2 flows and protostars that do not. In other 
words, the molecular outflow sources are no more - nor no less 
- clustered than the sources that don't drive H2 flows. This re- 
sult is confirmed by K-S testing: for the 6th NN radii, in the 
northern region there is a 71% probability that the protostar 
and H2-outflow-driving source radii are drawn from the same 
sample; in the south the probability is lower (11%) though still 
inconclusive - both distributions of NN radii are essentially the 
same. K-S tests for the 5th and 7th NN radii yield similar re- 
sults. 

We conclude that the differences we note at the beginning 
of this section are simply due to evolutionary effects. There 
are clearly regions within Orion A that are more evolved than 
others; young regions like L 1641-N that contain embedded 
protostars that drive H2 outflows, and older regions like L 1641- 
C with clusters of "protostars" that are not associated with large 
clouds, dense cores or molecular outflows. The high fraction 
of H2-outflow-driving protostars suggests that transition time 
scales - from protostar with H2 flow, via protostar without H2 
flow, to disk excess source - are relatively brief. The fact that 
the protostars without molecular outflows are no more widely 
distributed (on average) than protostars with outflows suggests 
that the protostellar phase as a whole is short in comparison 
to the overall time-scale over which young stars become more 
widely distributed. 

Finally, the NN analysis summarised in Table [TJ also 
demonstrates that there are considerable differences between 
OMC 2/3 and L 1641. All protostars, regardless of whether 
they drive H2 outflows or not, are far more clustered in the 
north (NN radii are consistently shorter) than in L 1641 in 
the south. Clustering, the distributions of protostars and disk 
excess sources, and the population of young stars through- 
out Orion, will be discussed in detail in a forthcoming paper 
(Megeath et al., in prep.). 
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Fig. 19. Histograms showing the distribution of 6th nearest 
neighbour radii for all Spitzer-identified protostars (open bars), 
for those that drive H2 outflows (black filled bars), and for those 
that do not drive H2 outflows (grey filled bars). Top - entire re- 
gion; middle - north of declination -6°; bottom - south of dec- 
lination -6°. 

5. Summary and Conclusions 

We present near-IR images of H2 outflows spread across 
an eight square degree region that encompasses most of the 
Orion A GMC. For much of this region Spitzer photometry is 
used to identify protostars and disk excess sources, and large- 
scale (sub)millimetre maps are used to plot the distribution of 
molecular gas. Tangential velocities (proper motions) are mea- 
sured for features in 33 outflows; these data are used to help 
associate the outflows with the protostars and dense molecular 
cores. 

We identify 43 new H2 outflows, increasing the number of 
H2 flows in Orion A to 116 (not including the Orion bullets). 



From the Spitzer sample of young stars we find sources for 
most of these outflows. Of the 300 or so protostars identified 
by Megeath et al. (in prep.), at least 1/3 seem to be driving H2 
flows. Indeed, the H2 outflow sources are predominantly proto- 
stars (Class 0/1 sources) rather than disk excess sources (mostly 
Class II T Tauri stars), though a few of the latter source sample 
do seem to power H2 flows. Most H2 outflow sources have pos- 
itive spectral indices and are associated with dust cores; most 
molecular cores that harbour protostars also drive molecular 
outflows. All evidence points to the extreme youth of H2 jet 
progenitors. 

We find no evidence for a preferred molecular outflow di- 
rection. Even in OMC 2/3, the region north of the ONC, where 
the molecular gas is confined to a north-south chain of dense 
cores, the H2 flows seem to be randomly orientated. Nor do 
we find a correlation between H2 flow length or opening angle 
and source spectral index or core flux at 1200yt/m. We find no 
evidence that "more evolved" sources, i.e. those with low a or 
weak 1200/im emission, drive longer or less-well-collimated 
H2 flows. The caveats associated with measuring flow parame- 
ters from H2 observations are discussed in light of these results. 

Finally, we investigate the spatial distribution of H2 out- 
flow sources in relation to the overall distribution of proto- 
stars. Although well-known star clusters, such as L 1641-N and 
L 1641-C, are associated with clusters of a dozen or so proto- 
stars and disk excess sources, only groups or clusters of pro- 
tostars that are still surrounded by dense molecular gas seem 
to be associated with multiple H2 flows. Protostars that power 
H2 outflows are no more (nor no less) clustered than protostars 
that do not. This, together with the fact that a high percentage 
of protostars do drive H2 flows, suggests that the H2-outflow- 
driving phase of a protostar's early evolution is only marginally 
shorter than the overall lifetime of the protostar; it seems clear 
that H2 jets fade very quickly as each source evolves towards 
the Class II pre-main-sequence phase. 
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Position of the H 2 outflow source if one is listed in column 4. Otherwise, the position of the brightest knot in the flow is given (enclosed in brackets). Note that the well-known sources of the 
HH 1/2, HH 83 and V380-Ori-NE flows were not identified in our tables of Spitzer protostars or disk excess sources. 

Most likely H2 outflow source from the Spitzer photometry together with the spectral index, a. Just a question mark means there is no obvious Hi flow source candidate. 
850 ^/m dust core coincident with the H 2 outflow source (from Nutte r & Ward-Thompson| 2007). The source must lie within a 14" radius (the JCMT beam at 850 yum) of the core position 
given by Nutter & Ward-Thompson in their Table Al. "Emission" or "undetected" means that no core appears in their table. However, "emission" means that the source identified in column 
4 is associated with diffuse 850/jm emission (surface brightness >100 mjy beam -1 , which is roughly equivalent to 5cr); a question mark means there are cores (or emission) in the vicinity 
which could be associated with the outflow or its source; "no obs" means that the outflow is outside the bounds of the SCUBA map. The core size (major/minor axis dimensions), Ap 850 , and 
integrated flux, Igso, are given for each core; Fgso is the flux, in Jy beam -1 , measured towards each IRS source (not necessarily the peak flux of the associated core). 

Same as for the 850 /im cores, but for cores identified from our analysis of the 1200 fim observations of Stanke et al. (in prep.). In this case the source must lie within 11" of the core centroid; 
"emission" corresponds to a surface brightness >75 mjy beam -1 (~ 5<r). The values given in brackets in the Ap l2 oo column are the orientation of the core major axis (this information was 
not available for the 850yum cores). 

Entire length (L) of all H 2 knots (in both lobes), or the distance from the source to the most distant Hi knot if only one flow lobe is identified; opening angle (6) measured from a cone centred 
on the outflow source that includes all H 2 features in the flow; position angle (PA) measured east of north. 
Associated HH objects, if any are known. 



Table 3. New H2 jets in Orion A (additions to the catalogue of Sta02). Columns are the same as for Table|2] 
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Appendix A: Description of the H 2 outflows in 
Orion A 

For a description of the new DFS outflows, and images of 
each flow, please visit http://www.jach.hawaii.edu/~cdavis/ 

Appendix B: Proper motion measurements 

For a description of the proper motion measurements and 
images of each flow with PM vectors superimposed, please 
visit http://www.jach.hawaii.edu/~cdavis/ 



This figure "davis-fg2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0812.3733vl 



